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Abstract 

Corresponding to optical Fresnel transformation characteristic of ray transfer matrix ele- 
ments (A, B,C, D), AD — BC = 1, there exists Fresnel operator F(A, B ,C, D) in quantum 
optics, we show that under the Fresnel transformation the pure position density \x) (x\ becomes 
the tomographic density \x} s ra (x\, which is just the Radon transform of the Wigner operator, 
i.e. 

f poo 

F \x) (x\ F f = \x) s rs>r (x\ = / dp'dx'S [x - (Dx' - Bp')] A (x',p') , 

J J — oo 

where s,r are the complex-value expression of (A, B,C, D). So the probability distribution for 
the Fresnel quadrature phase is the tomography (Radon transform of Wigner function), and the 
tomogram of a state is just the wave function of its Fresnel transformed state F^ \tp) , i.e. 
s , r (x\ ip) = (x\ F^ \ip) . Similarly, we find 

F \p) <p| F f = \p) 3 rs r (p\ = [ r dx'dp'S [p - (Ap' - Cx')] A (x',p') . 

J J — OO 

PACS numbers: 03.65.-w, 42.30.Kq 

1 Introduction 



In quantum optics theory, all possible linear combination of quadratures X and P of the oscillator 
field mode a and a* can be measured by the homodyne measurement just by varying the phase 
of the local oscillator. The average of the random outcomes of the measurement, at a given local 
oscillator phase, is connected with the marginal distribution of Wigner function, thus the homodyne 
measurement of light field permits the reconstruction of the Wigner function of a quantum system 
by varying the phase shift between two oscillators. In Ref. [1] Vogel and Risken pointed out that 
the probability distribution for the rotated quadrature phase Xg = [a^ exp(i9) + aexp(—i8)]/V2, 
[a, a^] = 1, which depends on only one 9 angle, can be expressed in terms of Wigner function, 
and that the reverse is also true (named as Vogel- Risken relation), i.e., one can obtain the Wigner 
distribution by tomographic inversion of a set of measured probability distributions, Pg (xg) , of the 
quadrature amplitude. Smithey, Beck and Raymer [HOH] also pointed out that once the distribution 
Pg (xg) are obtained, one can use the inverse Radon transformation familiar in tomographic imaging 
to obtain the Wigner distribution and density matrix. The Radon transform of the Wigner function 
is closely related to the expectation values or densities formed with the eigenstates to the rotated 
canonical observables. The field of problems of the reconstruction of the density operator from 
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such data is called quantum tomography. (Optical tomographic imaging techniques derive two- 
dimensional data from a three-dimensional object to obtain a slice image of the internal structure 
and thus have the ability to peer inside the object noninvasively.) The theoretical development in 
quantum tomography in the last decade has progressed in the direction of determining more physical 
relevant parameters of the density from tomographic data [51 El • 

In [51 [5] the Radon transform of Wigner function which depends on two continuous parameters is 
introduced. In this Letter we extend the rotated quadrature phase Xg to Xp = (s*a + ra^ + sa^ + r*aj /y/2, where 

| s | 2 — |r| 2 = 1, (s, r) are related to a classical ray transfer matrix^ M by 



x C D , 

s = i [A + D - i (B - C)} , r = -i [A - D + i (B + C)} , AD - BC = 1. (1) 

We shall show that the (D, B) related Radon transform of the Wigner operator A (x,p) is just the 
pure state density operator \x) s rs r (x\ formed with the eigenstates belonging to the quadrature 
xp. We name xp the Fresnel transformed canonical observable, or Fresnel transformed quadrature 
phase, so the probability distribution for the Fresnel quadrature phase is the Radon transform of 
Wigner function. These can be expressed neatly by 

F\x)(x\F* = \x) srsr (x\ = [ ( dx'dp'S[x-(Dx' -Bp')]A(x',p'), 



D = -( s + s*+r + r*), B = — (s* - s + r* - r) , (2) 

where Xp \x) s — x \x) s , and \x) s r — F \x) , \x) is the eigenvector of X e=0 , and F is the Fresnel 
operator corresponding to classical Fresnel transform in optical diffraction theory which we shall 
describe in the following. We name \x) s rs r (x\ the tomographic density operator. While the (^4, C) 
related Radon transform of A (x,p) is just the pure state density operator \p) s rs (p\ formed with 
the eigenstates belonging to the conjugate quadrature of Xp , 

F \p) (p\ Ft = \p) s rs r (p| = / ^ dx'dp'S [p - {Dx' - Bp')] A (x',p') , 

A = ~(s* -r* + s-r), C = i (s - r - s* + r*) . (3) 
2 2i 

Through ^ and (J3]) one can see how the quantum tomography is related to optical Fresnel transform. 

The optical diffraction transform is described by Fresnel integration whose parameters (A, B, C, D) 
are elements of a ray transfer matrix M describing optical systems, AD — BC = 1, M belongs to 
the unimodular symplectic group, the input light field / (x) and output light field g (x') are related 
to each other by Fresnel integration [lOl fTTj 



9 ( x> ) = ~nrm I 6XP 



— (Ax 2 - 2x'x + Dx'' 
2B v 



/ (z) dx. (4) 



In order to find the quantum correspondence (Fresnel operator) of Fresnel transform, the coher- 
ent state [HI [T3] \x,p) — exp[— -| (p 2 + x 2 ) + (x + ip) 0^/^/2] |0) is a good candidate for providing 
with classical phase-space description of quantum systems, thus we construct the following ket-bra 
projection operator 



i [A + D - i (B - C)] jj <inl P 



c *)(:)><(: 



(5) 



as the FO, where the factor ^JT. is attached for anticipating the unitarity of the operator F. In 
([5]) the symplectic transformation \ n n ) ( ^ ) mapping onto an Fresnel operator (FO) in 



Hilbert space is manifestly shown through the coherent state basis. Using the notation of \z) 
exp [— \\z\ 2 + za^~\ |0) , z — (x + ip) /y/2, and introducing complex numbers 



s = - [A + D - i (B - C)} , r = -- [A - D + i (B + C)] 



M a = i, 



(6) 



such that 



A B 
C D 



s —r 

-r* s* 



z 
z* 



exp 



-\sz — rz*\ 2 + (sz — rz*)a) 



|0>, 



so §5§ becomes [2] 



F(s,r) 



d 2 z 



\sz — rz ) (z\ 



Using the vacuum-state projector |0)(0| in normal ordering of boson operators 

|0) (0|=:e- Qta : , 



(7) 



(8) 



(9) 



and the technique of integration within an ordered product (IWOP) of operators [T^l HH [T7] we can 
directly perform the integration in ([5]) and obtain 



F{s,r) 



d 2 z 



ex P (' 



exp 



2s* 



Ii2ii2 I f\ *2 f 

|s| \z\ + sza' + z [a — ra'JH — — z H — —z — a' a 



exp 



— 1 a' a > : exp 



2s* 



(10) 



Note that this can be identified as a generalized squeezing operator with three real parameters 
[TBI HSI2D1- It then follows 



(z\F(s,r)\z') = 



1 



exp 



r ^ 9 r , z z 

z* 2 z 12 H 

2s* 2s* s* 



Then using the overlap between the coordinate eigenvector and the coherent state 

„2 



(x\ z) = it 1 ' 4 exp 



^ + V2xz-i^S 



(11) 



(12) 



and the completeness relation of coherent state we obtain the matrix element of F (s, r) in coordinate 
representation, 



(x'\ F (s, r) \x) = — (x'\ z) (z\ F, (s, r) / \z') (z'\ x) 



exp 



2B 



(Ax 2 - 2x'x + Dx' 2 ) 



(13) 



y/2mB 

which is just the kernel of generalized Fresnel transform in 

Now if we define g(x') = (x'\ g), f (x) — (x\ /) and using Eq. (|13[) . we can rewrite Fresnel 
transform in Eq. (JJ) as 



9) 



dx (x'\ F (A, B, C) \x) (x\ f) = (x'\ F (A, B, C) \f) 



(14) 



Therefore, the 1-dimensional FT in classical optics corresponds to the 1-mode FO F (A, J5, C) oper- 
ating on state vector |/) in Hilbert space, i.e. \g) = F (A, B,C) \ f). Using the Fock representation 
of \x) , 

lx 2 + V2xat- Vll°). (I 5 ) 



\ x) = tt 1 I A exp 



and flUl as well as (11211 we calculate 



d 2 z 

F (s, r) \x) = \fs I \sz — rz*) (z\ x) 

7T 



■ exp 



exp 



x 

~~2 



1, 




— \sz - 


- rz 


2 1 








z* 2 


M 2 


2 


2 



10} 



(16) 



where 



r-V4 



x/s* + r 



■. exp{- 



s — r x 



y/2x + S + 



7.2 



s* + r* 2 s* + r* 



s* + r* 2 



-}|o) 



(17) 



We can see that |x) sr — f(s,r) |x) is the eigenstate of Fresnel transformed quadrature phase, 
because from 

«l»>..r= f^-^T^Kr, (18) 



so 



s* + r* s* + r* 
[(s* + r*)a + (s + r)a t ] |i) ar = s/2x\x) ar . 



(19) 



This can be further confirmed by examining its completeness relation. Using the IWOP technique 
and © as well as n^prz — 2(r*+s+) + 2(V+s) ' we can P rovc tnat l x )ar m& ke up a complete set (named 
as the tomography representation), 



f 



dx \ x) 



dx 



+V2x 



\s + r\ s/it 



exp 



s — r s — r \ x x 

n> 



s* + r* s + r J 2 



a 1 a 



s + r a 1 s + r a 



s* + r* s + r J s* + r* 2 s + r 2 

Cfo; 1 / Q*r; -I- ?'nt -1- Q^;t 



s + r \/ir 



exp ■ 



s a + ra T + sa T + r a , , 
•'■ I J> : 



V2 



then using © and X = S ^-,P = as well as s* + r* = D + iB, s* - r* = A - iC, we can 

reform (|2H)l as 



1 = 



dx 



VD 2 + B 2 J-oo 



exp 



1 , ^2 

(x- DX + BP) 



D 2 + B 2 



and \x) is express as 



-1/4 



y/D + iB 



exp 



A-iCx 2 V2a; f D - iB a) 2 
D + iB~2 + D + iB a D + iB ~~2~ 



|0> 



and (1191) becomes 



(■D-^ — -B-P) \ x ) s ,r = x \ x ) 
According to the Weyl quantization scheme [21] 



H(X,P)= { I dpdxA(x,p)h(x,p) 

J J — CO 



where h (x,p) is the Weyl correspondence of H (X, P) , A (x,p) is the Wigner operator 

2 /"GO 



A(a?,i>) 



2tt 



X+ 2/\^2 



(21) 

(22) 
(23) 
(24) 

I, 

(25) 



we know that the classical Weyl correspondence of the projection operator \x) s rs (x\ is 



2nTr 



1 



— OO 

oo 



2nB ' dM6XP 
5[x- (Dx' - Bp')} 



ip'u + —u (x — Dx') 

B 



(26) 

From we see that under the Fresnel transformation l|16p the pure position density \x) (x\ becomes 
the tomographic density \x) s rs (x\, and further from (|24[) and ((26)) we see that it is just the Radon 
transform of the Wigner operator, i.e. 



Fix) (x\F^ 



dx'dp'5 [x - (Dx' - Bp')] A (x',p') . 



(27) 



Therefore, the probability distribution for the Fresnel quadrature phase is the Radon transform of 
Wigner function 

| (x\ Ft | 2 = | s , r (x\ ij}) | 2 = / r dx'dp'5 [x - [Dx' - Bp')] (^| A (x',p') |V) . (28) 



Moreover, the tomogram of quantum state \tp) is just the wave function of its Fresnel transformed 
state F' \ip) , i.e. s>r (x\ ip) = (x\ F' \ip) . (|27| and (p8)l are the new relation between quantum 
tomography and optical Fresnel transform, which may provide experimentalists to figure out new 
approach for generating tomography. 

Similarly, we find that for momentum density, 



F|p><p|F+ = b>,.„. r <p| 



dx'dp'5 [p - (Ap' - Cx')} A (x',p') 



where 



F\p) = \p) ( 



-1/4 



t2 



exp ■ 



D + iBp 2 y/2ip t A + iCa 
' A-iCY + A-iC a + A-iC~2 



|0). 



(29) 



(30) 



As an application of the relation (|27|) . we recall that the Fresnel operator F(r, s) makes up a loyal 
representation of the symplectic group, i.e. 



F(r,s)F'(r',s') 



d 2 z fd 2 z' 



: exp 



2s"* 



a* 2 I : exp 



1 ) a/ a > : exp 



2s" 



F(r",s"), (31) 



where 



or 



s r 
-r* s* 



Is"! 2 - lr"l 2 = 1, 



A" B" 
C" D" 



A B 



A' B' 



AA' + BC AB' + BD' 
A'C + CD B'C + DD' 



C D J \ C D' 
It then follows from ([27]) {33]) and that 

F'(r', s')F(r, s) \x) (x\ ft(r, s)F'V, s') 

dx'dp'5 [x- ((B'C + DD') x' - (AB' + BD') p')} &(x' ,p' 

i 

dx'dp'5 [x - (D"x' - B"p')] A (x',p') , 



(32) 
(33) 



(34) 



in this way the complicated Radon transform of tomography can be viewed as the sequential oper- 
ation of two Fresnel transforms. 

In sum, we have revealed the new relation connecting optical Fresnel transformation to Radon 
transformation in quantum tomography, i.e. the probability distribution for the Fresnel quadrature 
phase is the tomography (Radon transform of Wigner function). The tomography representation 
s ^ r (x\ is set up, based on which the tomogram of a state \ip) is just the wave function of its Fresnel 
transformed state , i.e. s>r (x\ tp) = (x\ . The group property of Fresnel operators help 

us to analyze complicated Radon transforms in terms of some sequential Fresnel transformations. 
The new relation may provide experimentalists to figure out new approach for realizing tomography. 
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